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ABSTRACT 

We present maps of the nature of single star progenitors of supernovae and 
their remnants in mass and metalhcity space. We find our results are similar 
to others but we have gone further in varying the amount of mixing and using 
various mass-loss schemes to see how the maps change. We find that extra- 
mixing, in the form of convective overshooting, moves boundaries such as the 
minimum mass for a supernova or WR star to lower masses. We also find 
that the pre-WR mass-loss determines the shape of our maps. We find that 
different mass-loss rates lead to quite different results. We find that the rise 
in luminosity at 2nd dredge-up places quite tight constraints on the masses 
of some progenitors and in particular the progenitor of supernova 2003gd. 

Key words: stars: evolution ~ stars: supernovae progenitors - stars: Wolf- 
Rayet 



1 INTRODUCTION. 



Supernovae (SNe) have most probably been observed since the dawn of astronomy. The 
earhest known surviving r ecord of a supernova was ma de in the eleventh century by ancient 



Chinese astronomers (see 



Stephenson fc GreenI (|2002| ) for the details of these early obser- 



vations), but the modern study of SNe began with Baade and Zwicky in the 1930's when 
they realised SNe are more luminous and rarer than the more common novae. The high 
luminosities and broad spectral lines led them to conclude that SNe are extreme explosions, 
the result of a dying star collapsing to a neutron star. This basic model still holds. Sixty- 



* E-mail: jje@ast.cam.ac.uk 



2 J. J. Eldridge and C. A. Tout 

five years of work lias vastly increased our understanding but many questions about SNe 
remain unanswered. 



A core-collapse SN occurs once a star has a core (usually of iron group elements) that 
cannot be supported by any further nuclear fusion reactions or electron degeneracy pressure. 
This leads to the core collapsing to a neutron star or black hole. In the process it releases a 
large flux of neutrinos that interact weakly with the surrounding stellar envelope depositing a 
large amount of energy (10^*^ J) that can drive the star to explode. In massive cores supported 
by electron degeneracy pressure core-collapse can be initiated by electron capture onto nuclei. 



Important questions are, 'Which stars give rise to which SNe and leave which type of 
compact remnant?' The answers will provide a test of stellar evolution theory, give informa- 
tion for galactic chemical evolution and predict SNe birth rates and remnant populations. 
To solve these questions we must use our understanding of stellar evolution to model stars 
up to the point at which they become SNe. Three main factors affect the answer, the initial 



mass and meta^ 
i ng its lifetime. 



licity of the star a nd mass-loss (or mass-gain in some binary systems) dur- 
Heger et al.l (|2003l ). using the models described in lWooslev. Heger fc Weaver 
( 20021 ) ■ have considered these factors and answered some of our questions. While their study 
is informative, it lacks an accurate metallicity scale with only three metallicities studied 
Z/Zq = 0, 10~^ and 1. They only investigate a single mass-loss scheme, do not investigate 
the lowest mass for a star to become a SN nor describe the SN progenitors. These factors 
limit comparison with observations. In this paper we expand upon their work and add de- 
tails on the progenitors of S Ne. These details ar e im portant for surveys of progenitors now 

Smartt et aJ (2003) and 



underway by various groups, 



Van Dvk et al. 



(20031). 



From the literature we see that there are a variety of mass-loss schemes in use. We have 
compared a number of these and present a useful parameterisation. Within our study we have 
found many regions of interest, Wolf-Rayet (WR) stars, red giants (RG) and asymptotic giant 
branch (AGB) stars. Their relative abundance, dependent on the initial mass function, and 
luminosity are of key interest as well as whether they do achieve the conditions to undergo 
core collapse. We have compared our results with observations of the known progenitors to 
date. In a few years there will hopefully be enough known progenitors to be able to fully 
test our predictions. 
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2 CONSTRUCTION OF MODELS. 

2.1 Stellar Evolution Code Details. 



We ha ve used the Cambridge STARS stellar evo l ution code originally developed bylEggleton 



f)l97l[ ) and updated most recently by 



Pols et al 



( 199,51) and lEldridge fc ToutI (|200|). Further 



details can be found at the code's home page (http:/ /www.ast.cam.ac.uk/~stars |. In this 
work we use 21 zero-age main-sequence models that have masses from 5 to 2OOM0 with a 
uniform composition determined by X = 0.75 — 2.5Z and Y = 0.25 + 1.5Z, where X is 
the mass fraction of hydrogen, Y that of helium and Z is the initial metallicity. It takes 
values from 10^^ to 0.05, equivalent to 2^00^® ^'^ The masses and metallicities used 

are shown on the axes of Fi gure 1. The models ar e evolved with a mixing length parameter 



a 



2.0 in the formalism of 



solar model (Pols et al 



Bohm-Vitensd (|1958[ ). This value gives a good fit to a realistic 



19981 ). 



We have also created maps with alternate compositions: X = 0.7, Y = 0.3 — Z and X = 
0.76— 3Z, Y = 0.24+2Z. With constant hydrogen abundance the variation with metallicity is 
accentuated. For example 2nd dredge-up occurs at higher masses with increasing metallicity 
than we find below. Changing the primordial hydrogen abundance from 0.75 to 0.76 as in 
the second case only produces small changes to our results. 



2.1.1 Nuclear Reactions and Opacities 
At the time of 



P"1s et al 



( 19951 ) only six composition variables were modelled. We have 



extended this to eight: H, He, C, N, O, Ne, Mg and Si. We have a lso updated the reaction 
rates according to the N ACRE collaboration ( Angulo et all 1999), replacing a number of 



those originally taken from 



Caughlan &: Fowled (|l988l ). We have further updated the reaction 



network to follow the later stages of burning although we encounter numerical difficulties 
at the onset of neon and oxygen burning because these reactions occur at a much more 
rapid pace than the previous evolution with timescales of a year and below. During these 
burning stages the envelope is only affected to a small degree since tnuciear << ^thermal- 
Therefore we use models after core carbon burning to estimate the details of progenitors. 
We have confirmed from models that do experience a neon burning phase that the late 
stages of evolution are rapid and have a negligible effect on the final observable state of the 
progenitor. 

We use the updated tables of lEldridge fc ToutI (|2004l ) that accurately follow the opacity 
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variation owing to the He/C/0 r atio changes during evo lution by including all the tables 



available from the OPAL project (jiglesias fc R,oger 



199(i ). We found important differences 



in the structure of Wolf-Rayet (WR) stars and so we are able to model these stars with 
increased accuracy. 



2.1.2 Convective Overshooting. 

Convective mixing is modelled simultaneously with the structure with a diffusion equation 
that allows us to correctly follow stable semi-convection. Our convective overshooting is not 
parameterised by a fixed fraction of the pressure scale height, Hp. Instead we modify the 
classical Schwarzschild criterion for instability, normally expressed as 

Vr > Va, (1) 

for convective instability to occur, where Vr and Va are the radiative and adiabatic gradients 
respectively. To include overshooting we alter this condition to 

V, > V„ - 5, (2) 
where 

~ 2.5 + 20C + 16(2' ^ ^ 

( is the ratio of radiation to gas pressure and 6ov is the overshooting parameter that we vary 
to include convective overshooting or to rer nove it. Further details of the inc lusion of this 
scheme in the Eggleton code can be found in lSchroder. Pols fc EggletonI ()l997l l. They found 
from fitting models to ecli psing binaries 6ny to be 0.12. A similar trend was confirmed with 



a large pool of binaries in 



Pols et al. 



( 19971 ). This leads to overshooting lengths /qv between 
0.25 and 0.32Hp in stars of mass range 2.5 to 6.5M0 for more massive stars the overshooting 
length remains around 0.3Hp. Another way of quantifying the effect of overshooting is to 
measure the mass of the convective core of a ZAMS model. For the 6.5Mq star the mass 
over which mixing occurs grows from 1.6 to 2.IM0, for a IOMq star the growth is from 3.2 
to 4.IM0 and for a lOOM© star from 83Mq to 87 Mq. 

In our calculations we use 5ov = 0.12 to include convective overshooting and set it 
to zero to remove overshooting. The amount of overshooting affects the core mass during 
evolution. This affects the nature of the remnants and also the mass of stars that undergo 
second-dredge up. The extra mixing owing to convective overshooting can be considered 
to mimic any mixing process such as rotation or gravity wave mixing, not just convective 
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overshooting. For a review of the various extra mixing process es and the possible anomahes 
of stellar evolution they could resolve see IPinsonneaultl ()l997l ). 



2.2 Mass-Loss Prescriptions. 



Mass-loss modifies the evolution of a star by affecting the surface conditions and the total 

mass of the star. While theoretical mass-loss rates do exist (Vink, de Koter & Lamers 2000, 

2001) it is sensible to use empirical mass- loss rates if we aim to model what is observed. 

We have tested the theoretical rates of Vink et al. and found the results agree well with 

those with empirical mass-loss rates. We divide our mass-loss rates into two main categories 

those for Wolf-Rayet (WR) stars, that have lost their hydrogen envelope, and pre-WR or 

hydrogen-rich evolution. The nature of mass loss changes once the hydrogen envelope has 

been removed and observations show quite different mass-loss rates. WR stars have greater 

mass loss than an OB star of the same luminosity. The mechanism for mass loss in OB stars 

is radiatively driven winds. The WR star mechanism is not kno wn for certain. It could be 

" II 

driven by radiation, pulsations or opacity ffHeger fc Langeiiil996) . 

Our resultant maps are strongly dependent on the mass-l o ss pre s cription. We have m ade 



a number of grids us ing prescriptions similar to Mevnet et all ( 1994 ). 



Heger et al. 



120031) and 



Drav fc Toud (|2003[ ). For ou r main maps we hav e adopted the prescription below, they are 
essentially the NL rates from lDrav &: Tout! (|2003l ) but with a modification at low metallicity 
10~^) as detailed below. In section 5.1 we present detailed comparison of the different 



mass-loss prescriptions we have tested. 



2. 2. 1 Pre- WR Evolution. 



We use the mass-loss rates of 



de Jager. Nieuwenhuijzen fc van der Huchtl ()l988l ). here after 



JNH, except for a small region at low metallicity and high stellar masses detai led below. 



The JNH rates are old and complex but considered to be the most accurate (|Crowthei 



20011 1 . The rates are dependent on surface luminosity and temperature and are derived from 



a large pool of observations. We apply the common 



M{Z) = M{Zq) X {Z/Zq^-^ as in 



Heger et al 



y ad opted s caling with metallicity of 



( 20031 ) and iDravl ()2003h . This scaling arises 



from the assumption that stellar winds are line driven and with lower surface opacity at 
lower metallicity there are weaker winds. However while there is agreement that mass-loss 
scales in this form there is a range of suggested values for the exponent. We have tested 
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the sensitivity to changes in this scahng by producing a grid with [Z/Zq)^'^ and higher 
values for the exponent (up to 0.69 suggested by Vink et ah As expected a larger exponent 
increases variation with metallicity: the type Il/Ibc SN boundary varies with metallicity to 
a greater degree. Previous authors have commonly assumed that this scaling applies at all 
stages of evolution even though there is some evidence that it varies with spectral type. This 
seems likely because the exponent assumes the winds are line driven while mass loss from 
giants is driven by some quite different unestablished mechanism. If we do not scale the 
mass loss for giant stars (Tefj < 3.7) the result is that the minimum initial mass for WR star 
formation decreases with decreasing metallicity. This is contrary to other results. In another 
test we used an exponent of 0.5 for all stars except hot OB stars (T > 10^) where we used 
an exponent of 0.7. This produces results very similar to those of our preferred results using 
0.5 for all stars. 



We encounter problems scaling the JNH rates when Z ^ 10 ^ for the highest mass stars 
(M ^ 15OM0 when Z = 10"^ and M ^ SOMq when Z = 10'^). The stellar parameters are 
outside of the region covered by JNH. Using simple extrapolation the mass-loss rate drops 
to very small values and leads to strange structure over the map. There is still a minimum 
mass for Ibc stars but some of the stars more massive than this minimum do not lose all 
their hydrogen before SNe. The problem can be seen in Figure ^ where with the JNH rates 
at low metallicity the minimum mass for Ibc SNe decreases with decreasing metallicity. This 
is because the nature of first dredge-up differs in these stars and complicates the mass- 



l oss h istorv. However, if we replace the JNH rates with those of iNieuwenhuijzen fc de Jager 
( 19901 ) (NJ) which are based on the same set of observational data but are expressed more 
simply a smoother structure is obtain. All the stars above the minimum Ibc SN mass do 
undergo type Ibc SNe. We also find that none of our models at Z = 10~^ undergo type 
Ibc SNe. We do not use the theoretical rates of Vink et al. because they do not extend 
down to these low metallicities. We consider the rates of NJ rather than JNH to be the best 
extrapolation method at the current time at low metallicity and high initi al mass. We have 



also experimented with the low metallicity theoretical mass-loss rates of 
and find similar results. 



Kudritzkil (|2002h 
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2.2.2 WR Evolution. 

WR evolution begins just before the hydrogen envelope is completely removed from the star. 
The stages of evolution run in the sequence of WN, WC and WO. In each case a particular 
element dominates in the emission spectrum, nitrogen, carbon or oxygen respectively. How- 
ever the_definition of when a model enters each stage is arbitrary. We use those defined in 
Drjjaild). 



• WN: When Xgurfacc < 0.4 and Teg > 10^ 

• WC: When X < 0.001 and 0.03 < {Xq + Xo)/Y < 1.0. 

• WO: When X < 0.001 and (Xc + Xo)/Y > 1.0. 



During WN and WC evolution we use the rates of iNugis fc Lamera (|2000| ) that depend 
on luminosity and surface abundance They find that the mass-loss rate depends strongly 
on composition as well as lumiriositv- This dependence gives lower mass-lo ss rates than the 
mass-dependent rates of iLangeii ()l989l ) and the mass-dependant rates that lNugis fc Lamers 
( 20001 ) also derive. We have tested both sets of mass- dependant rates and find that they lead 
to extreme mass loss from WR stars at and above solar metallicity and lead to the remnant 



Nugis fc Lamer j (120001) 



becoming a white dwarf. Our preferred rates are therefore the main 
rates. 

One straightforward test of this mass-loss prescript ion is to compare th e number ratios 



Drav fc ToutI 



p are tm 
I (20ol 



performed this 



of WR stars to O stars, and the WR subtype ratios, 
test in detail and found good agreement to observed ratios. The one exception was the 
WN/WR ratio with the theoretical value being too low. There are a number of plausible 
reasons for this such as binaries affecting mass-loss and misclassification of some WN stars 
due to changing spectral features as metallicity changes. 



3 SUPERNOVAE TYPE DETERMINATION. 

There are two main types of SNe those without hydrogen in their spectra, type I, and those 
with hydrogen in their spectra, type II. Type I progenitors have lost their hydrogen envelopes 
and are white dwarfs or WR stars, while type II progenitors have retained their hydrogen 
envelopes. This gives us a basic method to discriminate between them. 



8 J. J. Eldridge and C. A. Tout 
3.1 Type II 

For a star with a model at the end of core carbon burning we say it will become a type 
II SNe if there is any hydrogen in its envelope. Further more we know observationally that 
type IIP SNe have retained a good proportion of their original hydrogen. This hydrogen 
leads to the light curve following a plateau phase powered by a moving hydrogen ionization 
front. Type IIL SNe only retain a small fract i on of their hydrogen as the light curve decays 

( 20031 ) with type IIP SNe occurring until the 



linearly. We adopt the value of 



Heger et al 



mass of hydrogen in the star drops below 2.OM0 when type IIL SNe occur. 



3.2 Type I 



Type I SNe are divided into types la, lb and Ic. SNe la are extremely bright events thought 
to be the explosive carbon burning of a degenerate white dwarf that has reached the Chan- 
drasekhar mass (MciJ by accretion; this is not a core collapse event. However types lb and 
Ic are core-collapse SNe in WR stars that have lost their hydrogen envelopes in a wind or by 
binary interaction. They are discriminated by the presence or absence of helium lines in the 
spectrum. However the differentiation between models that are type lb or Ic is not obvious 
so we will label them as Ibc. All of our progenitors still retain some helium and the loss of 
this final layer is likely t o depend on the ve ry end phase of evolution and the explosion itself. 



We also adopt from 



Heger et al. 



( 20031 ) a method of determining the strength of a type 
Ibc SN. This will allow us to quantify the number of SNe that have no observable dis- 
play and are thus unseen. This allows us to estimate ratios of type II to type Ibc SNe. 
The largest progenitors are thought to give no display because the core is so massive that 
even with a large explosion energy nothing escapes the forming black hole. The exception 



woul d be if a iet driven SNe occured 



play (MacFadven. Wooslev fc Hegei 



:hat m akes a black hole and produc es an observable dis- 



20011 ). These ranges are taken from 



Ensman fc Wooslev 



(119881 ) and are similar to those used for determining the supernovae remnant (see next sec- 
tion). 

• Mcore(He) > I5M0, no display. 

• 15Mq > Mcore(He) > 8Mq, faint SN. 

• 8M0 > Mcore(He) > 5Mq, possibly faint SN. 

• Mcore(He) < SMq, bright SN. 
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4 SUPERNOVAE REMNANT DETERMINATION. 

Determination of the remnant formed at the heart of a SN is a black art. The physics is 
extreme and poorly understood. However there are many prescriptions which rely on the 
conjecture that a more massive core leads to a more massive remnant. We consider here 
three types of stellar remnant: white dwarfs, neutron stars and black holes. 

White dwarfs come in three flavours. They can be made of helium (a He WD), a mixture 
of carbon and oxygen (a COWD) or a mixture of oxygen and neon with some sodium and 
magnesium (an ONeWD). HeWDs are formed by very low mass stars or more probably in 
binary stars and are not of interest to us here. COWDs are also not formed in SNe but are 
the remnant cores of stars that lose their envelopes on the asymptotic giant branch (AGB). 
At solar metallicity stars in the mass range 0.8 to 7Mq (no convective overshooting) produce 
a COWD while larger stars, up to 9Mq develop an ONe core before the AGB. Some of these 
lose their envelope and form an ONeWD. The most massive undergo a SN driven by electron 
capture on to magnesium. The dividing line between these two is of great interest but is 
somewhat uncertain and depends on the mass-loss rate and the details of the thermal pulses 
in AGB stars. 

From our models we decide that stars with cores that grow by thermal pulses after 2nd 
dredge-up form COWDs or ONeWDs. This corresponds to a CO core mass < 1.38 after 2nd 
dredge-up. All other stars undergo a SN and form a neutron star or black hole. To distinguish 



between neutron stars and b 
compare that of 



ack h oles there number of different prescriptions. We 



Heger et al.l (j2nn3[ ) with ours 



4.1 He Core Mass. 



Heger et al.l (j2003|) use a simple system based on the helium core mass at the end stage of 
evolution. This is useful because the helium burning life-time is slower than for the later 
stages of burning so the core mass increases by only a small amount in the late burning 
stages. However the method relies on setting the values for different remnants by using 
simulations of SNe. The values used are: 

• Neutron Star: Mcore(He) < 8Mq. 

• Black Hole (by fall back): 8Mq < Mcorc(He) < ISM©. 

• Black Hole (directly): Mcore(He) > ISM©. 

These values are based on the SN models of iFrverl (|l999r i. By using this scheme we can 



10 J. J. Eldridge and C. A. Tout 



compare directly with iHeger et al.l (|2003l ). We should also note that SN that form a black 
hole directly also have no display. While those SN forming black holes by fall back will be 
faint. Unless a jet driven SNe occurs producing a black hole and an observable display. 



4.2 Core Mass by Energy. 

We use a core mass defined by the binding energy of the star. We assume a neutron star is 
formed at the centre of the star after core collapse of mass Mch = IAAMq. This produces 
about lO^^J of energy from the release of gravitational binding energy in neutron star for- 
mation. We then assume a hundredth of this energy is transferred into the envelope by some 
unknown mechanism. The current suggestion is the transfer occurs via neutrinos released 
from forming the proto-neutron star that are thermalised within the envelope or dense outer 
parts of the core. 

We integrate the binding energy of the star from the surface towards the centre until we 
reach 10^^ J. The envelope outside this region is ejected with the remaining amount forming 
the remnant. 

If we have a Mj-em > it is a black hole. The advantage of this method is that we 
are able to estimate the remnant masses and ejected masses quickly for large numbers of 
models. The resulting minimum mass for a black hole is similar to that from the He core 
mass method. However these remnant masses are at best a lower bound as we use a very 
simple model for a complex event. 



5 RESULTS. 



5.1 Comparing mass-loss schemes 

To demonstrate the effects of the main mass-loss prescriptions in use we evolve three stars 
at solar metallicity with six different mass-loss schemes. The prescriptions are detailed in 
Tabled 



The MM scheme is as used by iMevnet et all (jl99J) and HG is as used by iHeger et al. 
( 2OO3I ). The numerical constants are extra factors which are used to e nhance or supp ress 



the mass-loss in certain stage s of evolution. For example the rates of iLangerl ()l989| ) are 
reduced by | due to clumping (|Hamann fc Koesterkdll998h . For the HG WR rates we apply 
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Table 1. The mass-loss prescriptions identifiers. JNH and NJ are as described in the main text, NL is from lNugis fc Lamer j 
I2OOO) and L from Lanecr (1989J. 



Identifier 


Pro-WR 


WR 


NL 


JNH 


NL 


NN 


NJ 


NL 


iMM 


JNH 


L 


MM 


JNHx2 


L 


HG 


NJ 


Lxi 


3HG 


NJ 


L 



Table 2. Details of Models with different mass-loss prescriptions, Z = 0.02. Mi is the initial ZAMS mass, L is the star's 
luminosity and Mp the mass of the progenitor before SN. 



Mass- Loss 



M,/Mq 


Identifier 


i°gio(r^) 


SN Type 


Mp /Mq 


25 


NL 


5.53 


IIP 


16.0 


25 


NN 


5.40 


IIL 


9.73 


25 


iMM 


5.53 


IIP 


16.0 


25 


MM 


5.50 


IIP 


11.0 


25 


HG 


5.70 


IIL 


9.85 


25 


3HG 


5.70 


IIL 


9.85 



40 


NL 


5.36 


Ibc 


8.49 


40 


NN 


5.44 


Ibc 


9.43 


40 


iMM 


5.17 


Ibc 


5.41 


40 


MM 


5.13 


Ibc 


4.90 


40 


HG 


5.23 


Ibc 


7.02 


40 


3HG 


4.89 


Ibc 


3.25 


120 


NL 


5.31 


Ibc 


7.69 


120 


NN 


5.20 


Ibc 


6.53 


120 


iMM 


5.15 


Ibc 


5.11 


120 


MM 




none 


< 1.4 


120 


HG 


5.13 


Ibc 


5.23 


120 


3HG 




none 


< 1.4 



the scaling of IWellstein fc Langeii (|l999| ) with and without the reduction by one third. For 
pre-WR evolution JNH and NJ are based on the same set of observational data however the 
formulae given in the papers depend on 20 or 3 terms respectively. Thus JNH rates have 
much more structure over the HR diagram. For WR evolution NL rates are based on more 
recent data and depend on luminosity and surface composition while the L rates are older 
and depend only on mass of the star. 

From Table El it is possible to see many trends. Doubling the pre-WR mass-loss rate 
generally leads to smaller remn ants and less hydro gen i n the envelope of th e progenitor. Of 
special interest are the rates of iHeger et al.l (12003^ and lMevnet et all ()l994( ). In the sample 
above these schemes lead to smaller remnant masses and in the most extreme case a WD. 
For these models the evolution was halted once they reached the Chandrasakar mass. They 



occur because of high pre-WR mass-loss followed by mass dependent rates for WR evolution. 
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At lower masses the agreement with other prescriptions is better, ahhough the NN rates 
lead to a lower minimum initial mass to form a WR star. 

The overall shapes of our maps are consistent between mass-loss prescriptions but spe- 
cific values c hange between them. Of special interest is the comparison between NJ as in 



Heger et al.l 1)20031 ) and our preferred rates JNH. If we compare the two graphs in Figure Q 



we see that the former lead to greater mass loss and therefore a lower mass for Ibc SNe at 
solar metallicity. Our mini mum mass for type Ibc SNe with NJ rates is 4Mq lower at solar 



metallicity that that from 



Heger et al. 



( 20031 ) due to using the NL mass-loss rates for WR 
evolution. The NJ rates also produce greater variation with metallicity. This is because the 
NJ rates are extremely simple with only 3 terms in the formula for mass-loss while JNH 
employ 20 terms that can resolve greater structure over the HR diagram and produce more 
accurate models. The NJ scheme blurs out the details and makes the mass-loss too large at 
some points of evolution. Here we prefer to use the JNH mass-loss rates until at very low 
metallicities when the stellar parameters move outside the valid region of JNH where we use 
the NJ rates which can be extrapolated easily. 

From Table El it can be seen HG mass-loss rates lead to the production of lower mass 
progenitors. The result on our maps is to erode the structure in the upper right-hand corner. 
For example FigureElwith HG rates would have neutrons stars as the only possible remnants 
when Z ^ 0.03 and black holes only forming directly when Z ^ 0.004. 



5.2 Progenitor, SN and SN Remnant Diagrams 



Heger et al 



(20031). Although, owing to our 



Our maps tend to broadly agree with those of 
specific coverage in metallicity, we obtain more structure within the map. All of the structure 
is due to the mass-loss prescriptions we have used in producing these models. It is important 
to note however that while the broad nature of these maps are correct some of the fine 
structure would in practice be blurred out by different helium content, rotation, magnetic 



fields and other mixing processes. We assess our maps for each variable individually. 



5.2.1 Luminosity, Figure{^ 

Most important here is the 2nd dredge-up line on the left half of the diagram its position 
changes with inclusion of overshooting. Without overshooting 2nd dredge-up occurs for zero- 
age masses between 8 and 11 Mq. With overshooting this range drops to between 7 and 9 Mq. 
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Interestingly the most massive of these stars do form CO cores above Mch and go through 
carbon burning forming an ONe core. Then, before the surface convection zone penetrates 
the hehum core, a convection zone forms between the hydrogen and hehum burning shells 
that dredges carbon and oxygen into this intershell region thus decreasing the CO core 
mass. When the surface convection zone dredges up helium it also dredges up the carbon 
and oxygen. This behaviour therefore means that during 2nd dredge-up the surface carbon 
and oxygen abundances increase, rather than the normal decrease found in lower mass stars. 

Another feature in the diagram is the structure in the top right hand corner. As mass 
increases we see the last few red giant progenitors are quite bright but are then replaced 
by WR stars and the pre-SN luminosity drops abruptly. The luminosity increases again as 
higher mass WR stars form. At the highest initial masses these WR stars once more become 
less massive and therefore less luminous. This structure becomes less pronounced at lower 
metallicities owing to the lower mass-loss rates. 

5. 2. 2 Final & Ejected Mass, Figure 

The final mass diagram has very similar structure to the luminosity diagram and again 
reflects how mass loss is more severe in the high-metallicity region and drops off at low 
metallicity. The mass ejected by the SNe is also interesting because it demonstrates that a 
WR star ejects very little mass because it is more dense and so more tightly bound than a 
red giant progenitor. 

5.2.3 SNe Type, Figure]^ 

Again differences occur with the inclusion of extra mixing. The boundaries for the types 
of supernovae shift down by around 5 — 30Mq depending on the metallicity and boundary 
type. The type IIP/L boundary follows the type II/I boundary although the type IIL region 
becomes broader at low metallicity. The shape of the type II/I boundary is well defined and 
is similar between the two diagrams with no type Ibc and IIL supernovae occuring when 
Z = 10^^ in our range of initial masses. 

With Figures ID and El and the details of section 3.2 we can identify the mass regime 
where SN Ibc become unseen. Those SNe that form black holes directly will have no display. 
Unless a jet driven SNe forms a black hole and produces an observable display. Unseen SNe 
will affect the ratio of type Il/Ibc SNe observed. At low metallicities all type Ibc SNe are 
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unseen. We do not find bright type Ibc SNe in our sample as none of tlie WR stars liave 
lielium cores less than ^Mq. 



5.2.4 'S'-^e Remnant, figures\^ &\B 

In shape our maps are similar to those of 



Heger et al 



( 20031 ) although we do not go to 



zero metallicity. We also have more structure in the low-mass end of the diagram with the 
white dwarf to neutron star transition moving to a greater mass at higher metallicity. The 
change between the no-overshooting and overshooting case is minor with only the transitions 
moving down in mass while the shape remains the same. Noticeably at high metallicity and 
mass with overshooting it is possible to see that neut r on sta rs are formed rather than black 
holes. This is similar to the findings of Heger et al.l (2003 ) The transitions from neutron 

( 20031 ) although our value with 



Heger et a) 



star to black hole are at similar positions in 
overshooting is lower by 5Mq. 

Comparing to the energy method of determining the remnant we find the same transition 
point between neutron stars to black holes just below 25Mq without overshooting and 21Mq 
with overshooting. The structure is also quite similar although it is difficult to decide when 
direct black hole formation may occur rather than fall back on to a black hole. This would 
require a denser remnant and our IOMq remnant line seems to follow the He core mass 
remnant region for direct black holes fairly well. While looking at these diagrams it is 
se nsible to ask where different types of gamma-ray bursts (GRBs) may occur as discussed 
by iHeger et al.l ()2003r ). It is thought that GRBs require a progenitor with a small radius (a 
few Rq) so the jets that form from the material accreted on to the central black hole can 
punch their way through the stellar envelope. If the stellar envelope is too large the jets 
dissipate and lead to a normal SN. Therefore any of our models that loose their hydrogen 
envelopes to become WR stars and also form a black hole may also be the progenitors of 
GRBs. It will be interesting to determine how many stars that fit this requirement give 
rise to an unseen SN compared to those that give a bright display. This could lead to an 
observational check if it were known how many GRBs have an associated SN. 
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6 DISCUSSION 

6.1 Low-mass SN Progenitors (M < 25Mq) 

Stars with zero-age mass less than 25Mq give rise to red giant SN progenitors and should be 
the most common single star progenitors observed. The luminosity of these stars determines 
how many we can expect to find in current observational searches for progenitors. We find 
the luminosity of th ese obiects depends upon the occurrence of second dredge-up in the late 



stages of evolution (ISmartt et al.ll2002f ). Second dredge-up can increase the final luminosity 
from log(L/LQ) < 4.6 to log(L/L0) > 5.2. This increases the maximum distance at which 
the progenitors can be observed. However the details of second dredge-up are complex and 
depend on the amount of nuclear burning throughout the evolution of the star. The character 
of stars during second dredge-up depends on initial mass and metallicity and affects the final 
fate of the star and the material dredged to the surface. We are able to define five classes of 
red giant progenitors, the masses are for solar metallicity without convective overshooting. 

• M ^ 7Mq, stars undergo second dredge-up and thermal pulses with a central CO core 
as a thermally pulsing AGB star. These lose their envelope and leave CO white dwarfs before 
before their cores reach Mch- 

• M ^ 8Mq, during or after second dredge-up carbon ignition occurs in a shell because 
the core is degenerate and has a temperature inversion caused by neutrino losses. The star 
then undergoes thermal pulses with an ONe core. Since the core mass after dredge-up is less 
than Mch we assume these stars lose their envelope forming ONeWDs. These are Super- AGB 
stars. 

• M ^ 9Mq, carbon is ignited before second dredge- up (in a shell if the centre is degen- 
erate). Thus at dredge-up there is a growing ONe core. If this can reach Mch before the 
envelope is lost then the star undergoes a SN. The outcome will depend on the nature of 
the thermal pulses. 

• M ^ IOM0, the CO core is greater than Mch before dredge-up. However shell carbon 
burning (enhanced by a thin neon burning shell in the most massive stars of this type) drives 
a convection zone that reduces the size of the CO core to Mch so dredge-up can occur. This 
CO material is mixed with the envelope and increases the CO abundance at the surface 
during second dredge-up. After dredge-up the star has an ONe core of Mch which progresses 
to a SN by electron capture by Mg^^. 

• M ^ 11^0; the helium core or CO core masses are too great for dredge-up to occur. 
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The limiting mass for the hehum core is SMq and for the CO core I.SMq. Nuclear burning 
in these stars progresses until it cannot support the core and a SN occurs. 



The actual masses vary with metallicity and convective overshooting. As metallicity in- 
creases or decreases the transitions increase or decrease respectively by 2Mq at most. Convec- 
tive overshooting has a strong effect of lowering the masses by 3Mq. At solar metallicity the 
maximum mass for a star to undergo second dredge- up is SM©. With convective overshoot- 
ing the above five cases occur in the following mass ranges, M ^ SMq, GMq ^ M ^ 7M0, 
M ^ 7.5M0, M ^ 8M0 and M ^ QM© 

Another intriguing result of convective overshooting is that carbon burning ignites in 
smaller CO cores so, after second dredge-up, more of the stars have an ONe core smaller 
than Mch when thermal pulses begin and therefore produce more ONeWDs than without 
convective overshooting. 



Th ere are few sets of models to which we can compare ours. yRitossa. Garcia-Berro fc Iben 
( IQQol ). and references therein, detail similar models at solar metallicity for this range of 



Podsiadlowski et al. 



masse s. Compared to their models our mass limits are lower by a solar mass. 
( 20031 ) also detail stars that undergo second dredge- up and the implications for SNe and find 
similar limits higher that ours by a solar mass. This change is due to our updated opacity 
tables and nuclear reaction rates as our models with older tables produce limits that are 
higher by a solar mass. 



6.2 High-mass SN Progenitors (M > 25M0) 

The change between the types of SNe Il/Ibc is abrupt. However it is more difficult to 
differentiate between the sub types of Ibc or IIP/L. This is mainly due to the arbitrary nature 
of determining the way composition affects the SN dynamics. Our results have boundaries 
that vary less than previous authors because we apply the JNH mass-loss rates. We believe 
this is appropriate for single stars that are not rapidly rotating or in binaries. 

The problems encountered scaling mass-loss at low metallicity, Z ^ 10 ~^, mean that our 
SNe type boundaries become uncertain. Although the theoretical rates of Kudritzk] ( 2002 ) 
indicate that a change in the mass-loss rates in this region is expected. This agrees with the 
cut-off found in our maps. 
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6.3 Convolution with an IMF - SNe Ratios. 

Observations tell us that in SOa-Sb galaxies, the ratio of Type II supernovae over Tvpe Ibc su 



perno vae is between 1.2 and 12. In Sbc-Sd galaxies, it is between 2.1 and IG lCappellaro et al 



( 19971 ). At solar metallicity we find the ratio without overshooting to be 6.5, and with over- 
shooting to be 6.1 for a simple S altpeter IMF. This is slightly higher than the value of 



5 calculated by 



Heger et al 



(|2003l ). Generally the theory agrees with observation but this 
basic picture must be altered at lower metallicities because most type Ibc SNe become faint 
or have no display thereby raising this ratio. However binary systems are likely to dominate 
Ibc supernova production at low metallicity and this complicates prediction of the SNe type 
ratio. 



6.4 Remnant Initial Mass Population. 

Figure m shows the probability of a remnant having a mass Mrem- At the highest metallicities 
the remnants are all of very low-mass due to a maximum remnant mass cut off, with neutron 
stars dominant. Then as metallicity decreases the maximum remnant mass increases until 
the distribution becomes smooth with no maximum mass cut off. There is also a peak above 
the typical minimum mass for a black hole. These are the most common black holes at that 
metallicity. Although as the remnant mass used for these diagrams comes from a very simple 
calculation they are at best estimates of the distribution. 

Very massive (M ^ IOOMq) black holes can form in SN however this can only happen at 
low metallicity {Z < 10"'^). The population of such black holes will be very small as they are 
only formed in the most massive stars. At higher metallicities massive black holes are not 
formed as mass loss limits the size of the SN progenitor and thus the size of the remnant. 
Any black holes observed to be very massive in high metallicity environments must therefore 
have grown via accretion from a binary companion or must have formed from an an earlier 
population of stars. 



6.5 The Progenitor of 2003gd and other SNe. 



The progenitor of SN2003gd is the first red progenitor to be ob served. The observed lumi 



nosity is such that log^o(-^/-^o) = 4.3 ± 0.3 (jSmartt el al.ll20n4( ). This places the progenitor 



in the very faintest region of our maps. Figure |H1 shows the luminosity limits of 2003gd 
compared with the model luminosities at solar and half solar metallicity. The star cannot 
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Table 3. Comparing values from maps to observations of progenitors. Adapted from table in lSmartt et aJ 1200.^) . 
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have gone through second dredge-up. This places important constraints on the mass of the 
progenitor. However the luminosities of the candidate stars are both near to the upper limit 
of the observation so we cannot distinguish between the overshooting or non-overshooting 
case. The minim um m ass we can have for this star is SM©. 



Smartt et all (|2003l ) compiled details on the sample of well studied SN progenitors (shown 
in Table El and Figure IH)). We can calculate various details to check against these observa- 
tions. Most match well with our overshooting models although the various types of type II 
supernovae do not match our II P predictions. This may be due these stars being binaries 



for example recent observations (jMaund et al 



2nn4( ) have shown the progenitor to 1993J 



was part of a binary. T here are also strong arguments that 1987A was part of a binary 



too ( Podsiadlowski 



1992I ). Except for 1993 J the ejected masses also match up well to those 



calculated from observations. As this list grows in length and detail more tests against our 
maps will become possible. 



7 CONCLUSION 



We find that our results generally agree with lHeger et al.l ()2nn8l ). although our maps vary less 



with metallicity because we apply JNH mass-loss rates rather than NJ over most of the H-R 
diagram except for a small region at low metallicity and large initial mass. Observational 
evidence does match well with our maps and we can also point out that the most common 
type of progenitors, red giants, should be like 2003gd, faint compared to the rarer more 
massive stars. Therefore observations may find it difficult to find red progenitors unless 
some of those that undergo 2nd dredge-up do evolve to SN. In which case there will be a 
population of very luminous red progenitors that have not been observed to date. 
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Figure 1. Borders between SNe Type Il/Ibc. Dashed line for NJ rates without convective overshooting, solid line for JNH 
rates without convective overshooting. See Table^for description of rates. Generally the stars to the right of the lines are type 
Ibc SNe, those to the left are type II SNe. Although at Z ^ 0.001 for JNH rates stars to the right of the second line are type 
II SNe. The metallicities and masses or our models are used as labels on the axes. 




5 7 8 9 1011121J15 20^^ . ^'^ 100 120 150 200 5 7 8 9 101 11213 15 20^ 25 30 40 ^ 50 80 100^20 150 200 



Figure 2. Luminosity immediately before SN, left no-overshooting, right overshooting. Units on contours are log(L/LQ). The 
dot-dashed line represents the minimum mass for a SNe. 
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Figure 3. Left final mass of progenitor before SN, right mass ejected in SN. Both are with overshooting, units on contours are 
mass in Mq. 




Figure 4. SNe Type, left no-overshooting, right overshooting. Solid line - boundary between type II/I SNe, dotted line 
boundary between type IIP/IIL from total mass of hydrogen. The dashed-dotted line is the estimate of the location of the 
lowest progenitor mass for a SN. 
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Figure 5. He core mass, left no-overshooting, right overshooting. Dot-dashed line - inininmm mass for SNc and therefore the 
boundary between white dwarfs and neutron stars (Mcoro(He) = 1.4Mq), solid line boundary between neutron stars and black 
holes formed by fall back (Mcoro(He) = 8Mq), dashed line - boundary between black holes formed by fall back and directly 
(Mcore(He) = I5M0). Note this same diagram can be used to determine the visibility of Type Ibc SNe. There are no WR 
progenitors of M ^ SMq, but the solid line separates possibly faint SNe and faint SNe, while the dashed line separates faint 
SNe from SNe with no display. The SNe forming black holes directly are those that have no display. The faint dotted line plots 
the boundary of type II/I SNe so the number of SNe with no display can be estimated. 




Figure 6. Remnant mass by Energy method, left no-overshooting, right overshooting. Units on contours are mass in Mq. 
Dot-dashed line represents the minimum mass for SNe. 
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z = 0.02 z = 0.01 




Remnont Moss / M© Remnont Moss / Me 



Figure 7. Graphs of probability of remnant masses obtained by convolution of the remnant mass diagram with overshooting 
and a Salpeter IMF. The value at the top of each individual graph is the metallicity plotted. The probability has been divided by 
the mass of remnant. There are sharp cut-offs at lower metallicity and the distributions become smoother at lower metallicity. 
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Figure 8. Comparing observations of 2003gd with models at solar and half solar metallicity. Both cased with or without 
overshooting fit the observations however the mass range is tightly constrained in both cases. 




Figure 9. Comparing values from maps to observations of progenitors. The contours are as in Figure 2 while the thick lines 
arc from F igure 4. The map is with overshooting and zoomed in on the region where observations exist. The data points are 
taken from ISmartt et . 



3 J2003t) . Filled circles represent type IIP SNe, the circle with filled centre represents the type IIL SN, 
the triangle represents the type Iln SN, the inverted triangle represents the type lib SN, the star symbol represents the Il-pec 
SN and the square represents the type Ic SN. 



